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Abstract 



Self-dual solitons of Chern-Simons Higgs theory are examined in curved 
spacetime. We derive duality transformation of the Einstein Chern-Simons Higgs 
theory within path integral formalism and study various aspects of dual for- 
mulation including derivation of Bogomolnyi type bound. We find all possible 
rotationally-symmetric soliton configurations carrying magnetic flux and angular 
momentum when underlying spatial manifolds of these objects comprise a cone, 
a cylinder, and a two sphere. 
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1 Introduction 



Since the Chern-Simons gauge field made appearance in physics literature as a mass term in 
three dimensional gauge theory [I]], it has attracted the attention in various ways, e.g., anyon 
and fractional statistics 0, a relation to knot theory [||, and conformal field theories in two 
dimensions Q . Chern-Simons gauge theories coupled to scalar matter field also support the 
soliton solutions carrying fractional angular momentum [[5], |6|. An example which provides 
static multi-soliton solutions is relativistic Chern-Simons Higgs model of a specific <j) 6 scalar 
potential, and it admits both topological and nontopological solitons |7|, [8], [| . For the vortices 



in Abelian Higgs model with critical coupling, it was proved in Ref. |10[ that the Bogomolnyi 
type bound is saturated after the inclusion of Einstein gravity. The progress to this direction 
has also been achieved in Chern-Simons Higgs model coupled to background gravity [O and 
Einstein gravity |fL2|| . The solitons of Abelian Higgs model are distinguished from those of 
Chern-Simons gauge theories by the fact that whether they are spinning objects or not. 
This angular momentum affects the derivation of Bogomolnyi bound in curved space: the 
Bogomolnyi limit of Nielsen-Olesen vortices is obtained under the static metric and the same 
4 scalar potential as in flat space case, however that of the self-dual Chern-Simons solitons 
is saturated under the general stationary metric and thereby introduce a <fr s scalar potential 



of negative coefficient in addition to the cf> potential of flat spacetime [0 . Furthermore, the 
scalar potential in the Bogomolnyi limit of the Einstein Chern-Simons Higgs model includes, 
at least, one parameter, even if we use the condition that the obtained first-order equations 
reproduce original second-order Euler-Lagrange equations and their solutions carry finite 
magnetic flux (or equivalently charge) and angular momentum. 

In this paper, for the static configurations of Chern-Simons Higgs model, we shall derive 
conventional Bogomolnyi type bound of the total energy defined by Euler number in terms 
of magnetic flux, the topological charge of the system. Examining the obtained Bogomolnyi 
equations by use of both analytic and numerical methods, we shall show that there exist 
only two types of smooth solutions despite the complicated 8 potential unbounded below. 
They are topological vortices and nontopological solitons which coincide exactly with the 
solitonic spectra in flat spacetime. First, when boundary value of the scalar field has a Higgs 
vacuum, self-dual vortex solutions are supported and the underlying manifolds of these 
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solutions are asymptotic cone and cylinder. Second, when boundary value of the scalar 
field has a symmetric local minimum, self-dual nontopological solitons are produced and the 
corresponding spatial manifold constitutes two sphere, asymptotic cylinder in addition to 
asymptotic cone. It is an intriguing point that the scalar potential vanishes at the boundary 
value of both solitons in open spatial manifold, which implies zero cosmological constant at 
asymptotic open space. Examining asymptotic behaviors of solutions for multi-solitons, we 
conclude that there exists an upperbound for the vorticity to form a vortex and the obtained 
solutions are decaying fast as the radial distance increases. 

A way to envisage the role of topological excitations in the path integral formalism of 
Abelian gauge theories is to reformulate the given theory through the dual formulation. 



It has been studied for the lattice version of the Abelian Higgs model fl~3]|, and for the 



continuum version of the Chern-Simons Higgs model ||14|| . In this paper, we shall rewrite 
the Chern-Simons Higgs model coupled to gravity in terms of dual gauge field. In the dual 
transformed theory, we obtain the explicit forms of the nonpolynomial interaction between 
the dual gauge and Higgs fields, and the topological interaction between the dual gauge field 
and the topological sector of the scalar phase. Furthermore, we demonstrate a role of those 
interactions in order to produce nonperturbative excitations and their mutual interactions. 
For gravitational field, the path integral measure of the duality transformed theory possesses 
a Jacobian relative to that of original theory. 

In the next section, we will introduce the Chern-Simons Higgs model in curved spacetime 
and rederive the Bogomolnyi type bound. In section 3, we rewrite our model through the 
duality transformation and discuss physics of the topological excitations comparing with the 
original theory. In section 4, we examine the Bogomolnyi equations and obtain all possible 
regular rotationally-symmetric solutions. We show that the decaying property of solutions of 
the Bogomolnyi equations and nonexistence of solutions with prescribed vortices in section 
5. Conclusions with some comments about our results are presented in section 6. 



2 Self-dual Chern-Simons Solitons in Curved Space 

In this section, we recapitulate the derivation of the Bogomolnyi bound for Chern-Simons 
solitons in curved space. In order to obtain such Bogomolnyi-type bound, we consider the 
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energy defined by the spatial integration of two dimensional scalar curvature instead of the 



matter action for static configurations as has been done in Ref. [12 . 

The action for an Abelian Chern-Simons gauge field theory with Higgs mechanism is 
written in curved spacetime 



S = J d 3 x-y/g 



1 „ « w . „ , 1 



R + -z—pA&At + -g^D^D^ - V(\<j>\) 



(2.1) 



16nG 2 " " " T 

where = |0|e if \ gauge covariant derivative is = — ieA^, and scalar potential V(|0|) 
will be fixed to give a suitable Bogomolnyi type bound. Dynamics of gravitation is governed 
by 2+1 dimensional Einstein gravity and the metric takes the stationary form, which is 
compatible with static spinning objects 

ds 2 = N 2 (dt + Kidx 1 ) 2 - j ij dx i dx j . (2.2) 

In the above equation the components of the metric, N,Ki, and 7^, = 1,2), do not 
depend on time variable. 

Variation of the action with respect to the fields leads to the equations of motion 

-±=_D, {Jgg^DA) = (2.3) 



k e 



fiup 



2 ^ P -ef, (2.4) 

G, lu = 8nGT^. (2.5) 
In Eq. (|2.5| ) denotes symmetric energy-momentum tensor 

Tpu = ^(D^D^ + D^D^-g^^g^D^D^-Vy (2.6) 

This theory possesses a local U(l) gauge symmetry so that it admits a conserved U(l) charge 

Q = Jd 2 x^f } (2.7) 

where 7 = det7ij and j is time- component of the U(l) current 

f = - % -^D^-D^^). (2.8) 

We are interested in the static soliton-like excitations satisfying the classical equations of 
motion ( p.3|) -( ^5|) , and particularly the self-dual solitons supported by first-order Bogomolnyi 



equations. Though it is in general difficult to have the expression for total energy of the 



system, we take it as an Euler invariant [19, 20]. From time-component of the Einstein 
equations (|2.5|), we have the energy expression in terms of matter fields 



E 
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\cl d2x ^ 2R 

d^l^AM 2 + ^D^Drf + V{\<j>\) - 



-N 2 K 2 



32nG 



(2.9) 



where 2 R is two dimensional scalar curvature, Di<p = (di—ieAi+ieKiAo)<j), K = diKj, and 

e*- 7 is two dimensional Levi-Civita tensor density of e 12 = ei2 = 1. Suppose that the potential 

I 

V is made of both G-independent U and G-dependent W such as V = —U 2 + W . We can 



rearrange the terms in the integrand of right-hand side of Eq. 
N{x % ) = 1 as follows 

U ^ 2 



under an assumption 



E = Jd 2 x^l [ j\<P\ 2 (A ± 



1 



+ - 7 iJ (A0 T i^e ikl kl D l( f))(D j( j) t i^e jml mn D n ^ 



T 2 



-v 2 ) 



TeA \ 



—d % {A 3 - K 3 A ) + KAq + -A \ 

,2 



U 



+ 



2 v v z 



2k 

2 ' A 



2 -v 2 ) 



(2.10) 



ev 



32nG 



± $ 

2 



+ / d 2 x <9, 



2 1 
■ ■ / PI) PD 1 — ~ 

& {±^KMo =F 6 —C) ± i (0 J D J - </>D. 



where $ is the magnetic flux defined by 

$ = -~ J d 2 x e<%, 



(2.11) 



and C a constant introduced so as to make the surface term on the last line vanish. Since 
the spinning objects specified by the nonzero angular momentum J given as 

1 



J 



8nG J\x\ 



dx l Ki 



(2.12) 
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are of our interest, the requirement that the surface term in the last line of Eq. (|2.10|) has 
no contribution to the energy relates C to the value of scalar magnitude at spatial infinity 
0oo = lim — such as 

>oo V 

C =!-<&,. (2.13) 



The offdiagonal metric components K in Eq. ( |2.10|) are replaced by matter fields from Oi- 
components of the Einstein equations ( f2.5| ) 



(A 2 - s^l D\ 

K = -SttG^ 1 N f \ (2.14) 

where D is an integration constant. 

Now let us choose a scalar potential to achieve a Bogomolnyi-type bound. Specifically, U 
is chosen to make the square bracket on the fourth line of Eq. (|2.10|) zero, and thereby the 
condition, that the first line of Eq. ( |2.10| ) should vanish, forces the auxiliary field Aq to be 

A = T -^ = T^(\<i>\ 2 -v 2 ). (2.15) 
e\d>\ 2k 



The third line of Eq. ( [2.10|) does not contribute to the energy because of the Gauss' law 



given by time-component of the Chern-Simons equations (|2.4j ). We read G-dependent terms 
of the scalar potential W as an eighth-order term from the fifth line of Eq. (|2.10p . Hence, 
for the configurations to satisfy the first-order equation 

A0 =F i^e ij ^ jk D k 4> = 0, (2.16) 

the energy is proportional to the magnetic flux and the Bogomolnyi-type bound is saturated 
if such solitons carry magnetic flux. The Gauss' law says that these flux-carrying solitons 
are charged objects and nontopological solitons can also be produced when the meson mass 
in the symmetric phase of the theory is not larger than e 2 v 2 /2n. The remaining equations 
are r/'-components of the Einstein equations ( |2.5| ) summarized as follows 

1 



87rGiV 



- V J V J )iV 



2 im \N ed>J\N e\d>\) J \ 32nG 



+ 1 



1 \H 



+ 



+ {D^ ± i — nfoDtyifiV =F l — lmnD" 

{D*(f> ± i^ lkl &4>){&<j) T f—lmnD^ 



a/7 x/7 



(2.17) 

Note that any self-dual soliton solution of the Bogomolnyi equations with iV = 1 automati- 
cally satisfies the above Einstein equations. 

From now on we fix the conformal gauge for 7^ 



and the Coulomb gauge for K\ ViK % = 0, which makes K % be expressed as follows 



K* 



k e 



v 



-_dj ln-0. 



(2.18) 



(2.19) 



e 2 v 2 3 

Substituting Eqs. ( |2.18|) -( p.l9| ) into 00-component of the Einstein equations ( p.5|) and re- 
placing the gauge field to the scalar field by use of the Bogomolnyi equation ( |2.16| ) 



A % - KiA = -(ditt =f e ij d j In 
e 

we obtain an expression of the metric function b as 



(2.20) 



^h{z)+h{z 



G 



rip=i \ z z p 



? 12 



(2.21) 



In the above h(z) {h(z)) is a holomorphic (an anti-holomorphic) function and the variables 
with tilde are dimensionless quantities 



z = x 1 + ix 2 



e 2 v 2 



(x l + ix 2 ), f 



G = AnGv 2 



(2.22) 



\K\ V 

If we eliminate the gauge field and the metric b by use of Eqs. ( 2.16 )- (|2.21 ), we have a 
Bogomolnyi equation from the Gauss' law which is the time-component of the Chern-Simons 
equations (|2.4j) for the gauge field 



d 2 \nf 



Ilp=i \ z z p\ 2 



G 



f 2 -j{f-2f 2 + 



(f ~ 1) 
=F 2e ij didjn, 



(2.23) 



and the equation for if) from Eq. ([2.14 ) 



= -^W) (^^Z^) G (f 4 - 2 f 2 + > (2.24) 
* \ llp=l \ z z v\ ) 

where d 2 is flat-space Laplacian, i.e., we will use d 2 = = A. 



The undetermined constant D in Eq. (|2.14 ) is fixed by the requirement that the above 



equations (|2.23|) - (|2.24|) reproduce the scalar equation (|2~3| ), D = 1 — (j) 2 ^. In synthesis, the 
scalar potential becomes 

Though the Bogomolnyi limit of a theory selects usually a unique scalar potential in curved 
spacetime fllQl , shape of the scalar potential ( |2.25| ) depends on the boundary value of the 



scalar field 0oo (see also Fig. [T]). This potential includes G-dependent eighth-order potential, 
and it is not bounded below since the coefficient of the eighth-order term is negative for 
positive Newton's constant G as shown in Fig. [I]. However the energy in order to support 
flux-carrying solitons of the Bogomolnyi equations is positive definite. Inserting the metric 
( |2.21| ) into the Euler characteristic ( |2.9| ), we have 



AkGv 



l —\ [d 2 xd 2 \nf\\z-z p \ 2 - [d 2 xd 2 \n\4>\ 2 + [d 2 xd 
Gv 2 [J p x =1 J J 



AttGv 2 \J ™~ „ X = V~ ~ v ' J ™~ ' J ~~ " v 2 

Jd 2 xd 2 \nip\. (2.26) 



From the above expression there can exist contributions from the second and fourth terms 
besides the first term for n ^ if the Bogomolnyi equations ( 2.2cl )-( p.24j ) contain the finite 



energy solution which behaves as \<f>\ ~ |x l | _e (e > 0) for large \x\, and ip ~ |x J | _2a (a is an 
arbitrary number) when ^ 1. We shall present the detailed analysis for the existence of 
such solutions in section 4. 

Looking at shapes of the scalar potential in Fig. [I], one may raise an intriguing question 
in connection with local supersymmetry in 2+1 dimensions. According to Ref. [p^] , a well- 



known argument in 2+1 dimensions is that local supersymmetry can ensure the vanishing 
of the cosmological constant without requiring the equality of boson and fermion masses. 
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Figure 1: Schematic shapes of the scalar potential V(|0|) for e 4 /8/t 2 = 1 and irGv 2 = 1. The 
solid, dashed, and dotted lines correspond to 0^ = 1, 0, and Jl + 1/27^2. 
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Suppose that the above theory is bosonic part of a presumed supergravity theory as has 
been done for N = 2 supersymmetric theory of self-dual Chern-Simons Higgs model in flat 



spacetime ||21|| . It is questionable that whether or not the supergravity version of our model 
satisfies the above criterion because of the following reasons. Although we let the cosmo- 
logical constant vanish from the beginning, the scalar potential ( |2.25| ) is unbounded below 
and has negative local minima. Specifically, degeneracy between the symmetric vacuum, 
|0| = 0, and the broken vacua, \<f>\ = v, disappears because of G-dependent terms, and this 
phenomenon is guaranteed by global supersymmetry [21j (see Fig. |IJ). In more detail, the 
value of the scalar potential for the boundary values of soliton configurations are calculated 
as follows; 

VM = l^oo) = l^T-fctfao - I) 2 (Voc - ^) • (2-27) 



If the boundary value of the scalar field is different from 0, 1, 1/vnGv 2 , then ^(0oo) 7^ 
and we may suspect that our model is not embedded in the 2+1- dimensional supergravity 
models claimed in Ref. [2^]. In section 4, we will show that all possible regular self-dual 
solitons with rotational symmetry have boundary conditions (f)^ = or 1 so that the model 
of our consideration shares properties of the vanishing cosmological constant, i.e., V(\<f>\ = 
O;0oo = 0) = V(\<f>\ = v ; 0oo = 1) = 0. One more comment should be placed: For both 
cases to support regular self-dual solitons (nonperturbative spectra), the potential change 
due to G-dependent terms does not affect the mass spectra of perturbative gauge and scalar 
bosons. Therefore, equality of the gauge boson mass and the Higgs mass in broken phase 
(0oo = 1), still holds in curved spacetime, i.e., m gauge = mni ggs = e 2 v 2 /\n\. The above 
properties imply that one of our present systems classified by an undetermined parameter 
0oo may be obtained as the bosonic sector of some supergravity models. 



3 Bogomolnyi Bound in Dual Formulation 

We will construct duality transformed theory of our interest and see how the solitons arise 
in this framework. Let us work in the path integral formalism 

Z = Jidg^dA^Wmd^V 3 , (3-1) 
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where the action S was defined in Eq. Q2.1| ). The interaction between the scalar and gauge 
fields is linearized by introducing an auxiliary field as follows 

exp J Sx^-g^^n-eA^n-eA^ 

= J [*3J II ^ ex P < / rf3 * {-^p^ + sTCitfJl - eA u ) | . (3.2) 

An auxiliary vector field is classically identified with the U(l) current by constraint 
equation. Contribution of the scalar phase Q(x) is divided into two: one from topological 
excitations described by a multi- valued function Q(x) and the other from fluctuations around 
superselected topological sector expressed by a single- valued function r)(x). Then the path 
integral measure for the scalar phase becomes 

[d£l] = [de][drj\, (3.3) 

and 77-integration explains conservation of the U(l) current in the path integral formalism 

J [dr,] exp J d 3 x y/gfTC^ri} « -L^C). (3.4) 



Dual vector field is introduced such as 



J [dC^SiVpC) ■■■ = J [dH^dC^i^gC - -^d v H p ) ■■■ , (3.5) 

where the scale n/e is introduced for later convenience. The action at this stage is quadratic 
in both and so that the path integrals are Gaussian type. Integrating both A^ and 
C M in a closed form, we obtain duality-transformed theory; 



igtdg^idH.m-'dmd®} 

r ( 1 n 2 k e^ iup 



K f pyp 1 



+-.— H ?»OpQ + - V(\<f>\) , (3.6) 



where = d^H v - d u H^. 
Equations of motion read 



1 dV k 2 

^(vurwn = - M + ^m^, (3.7) 
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= 8nGT°, (3.9) 

where the energy-momentum tensor is 



K 

e 



2 



T° = d,\(f>\dM-^-r 2 9 p(7 H, p H ua (3.10) 



In the Higgs phase with \<f>\ = v, the interaction term between the Higgs field and the 
dual gauge field turns out to be Maxwell term so that the Chern-Simons Higgs theory is 
transformed to topologically massive gauge theory and thereby both theories support an 
odd-parity helicity-one photon of mass e 2 v 2 /\n\ [|1^]. Even though the action in the duality- 
transformed theory ( |3.6|) looks inappropriate to describe the symmetric phase of Chern- 
Simons scalar electrodynamics because of the singular factor 1/|0| 2 in both the action and 
the path integral measure, this term plays an important role to produce nonperturbative 
spectra, e.g., nontopological solitons. For this, let us look at the Gauss' law in flat spacetime 



*(i!K B - ± ¥|r<* , -4>- « 3 - u) 



Since the electric field term in the left-hand side of Eq. (p. 11 ) has no long-range contribution 
in the Higgs phase, the magnetic flux must be quantized. If there exists a soliton of which 
the scalar field behaves like \<ft\ ~ \x l \~ e (e > 0) for large \x l \ in addition to the discrete 
magnetic flux, the magnetic flux has additional contribution to compete with that of the 
electric field term (|3.11|) , which is not necessarily discrete. It is indeed the case of the Chern- 
Simons Higgs theory |5) . Since the gauge coupling e is inversely coupled to the Maxwell-like 
term, which describes interaction between the scalar and gauge fields, the strong coupling 
expansion is achievable when the nonpolynomial Higgs interaction is neg lected |]|. Note 
that, though the classical gravity is not affected by the duality transformation, path integral 
measure for the gravitational field contains a nontrivial Jacobian factor which comes from 
gauge dynamics |T7| . 

From now on we explore the static solitons given by solutions of the classical equations of 
motion. The duality transformation makes the system complicate by changing the first-order 
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Chern-Simons equation ( ^Jj) to the second-order equation ( p78|) , however Eq. ( p78|) is reduced 
to the first-order Chern-Simons equation if we replace the gauge field by the dual gauge 
field H^. It is due to the fact that a relation between the auxiliary field and the dual gauge 



field H^ in the right-hand side of Eq. (|3.5| ) is the form of classical Chern-Simons equation. 
Instead of solving the Euler-Lagrange equations (|3.7|) -(|3~9|), we derive a Bogomolnyi-type 
bound of the duality-transformed Einstein Chern-Simons Higgs theory again. Similar to the 
previous section, we define the energy as Euler invariant and assume the 00-component of 
stationary metric N 2 to be 1. Reshuffle of the terms of the energy gives us an expression of 
the energy such as 

( ,^2 



E 



D 



2e 2 |0 


I 2 




K 


4( 


e\</>\ 


e 2 


e ij 


' 2k 




K 

+ — 


H 



Y J dAH ± — 

v 2k 



V 



H, 



e e 



'.J 



ij 



„2 



r C 



d i (K J (-H + 



H 



(3.12) 



ev 2 



+ (w + -k(± 

d 2 x e ij H 4 i 



v 2 )H n 



+ / d 2 x di 



ev 



ev 



±—<»K j (H* =F -^Cn)* =F -^d 5 H Q 



4 

ev 

— <^D 

Ik 



K 



32nG 



2e 



rJJ 



-.{Hij + KidjH — KjdiH Q ), and Co is an integration constant. 



where H — ••,<— v-<-<-m <.»,^,<.<u <>,■ 

2^7 J 2^7 V J J u 3 

Similar to the procedure done in the section 2, we obtain the Bogomolnyi-type bound for 
the duality-transformed Einstein Chern-Simons Higgs theory. Square brackets in the second 
and third lines in Eq. ( |3.12| ) vanish by use of the equations for the gauge fields (|3T8|). Let us 
assume that configurations satisfy 



H = =F — 
^2k 



(3.13) 



Then we can fix Co as 1 — 0^ from the condition that the objects of our interest have spin 
J and the boundary terms in the last line of Eq. ( |3.12| ) vanish. With the aid of the solutions 
of Oi-components of the Einstein equations, we choose a specific form of the scalar potential 
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to let the terms in the fourth and fifth lines be zero. Hence the Bogomolnyi-type bound is 
attained for the solutions of the following equation 



H = ±^U = ±%- I0IW-A (3-14) 

Replacing the gauge and gravitational fields in Eq. (|3.14j) by use of the relations obtained 
above, we have the same Bogomolnyi equations in Eq. ( |2.23| ) and Eq. ( [2.24 ). Among the 
remaining equations, solutions of the Bogomolnyi equation satisfy r/'-components of the 
Einstein equations (|3.9| ) and reproducibility of the scalar equation (|3.7|) determines the in- 
tegration constant introduced by Oi-component of the Einstein equations. 



4 Soliton Solutions 



In this section we explore soliton solutions by examining the equations of motion, i.e., one 
is Bogomolnyi equation fl2.23|) and the other is the remaining Einstein equation ( 2.24 ). We 
restrict our interest to rotationally symmetric configurations satisfying the following three 
conditions; 

1. nonsingular solutions of equations of motion, 

2. solutions of finite energy which is proportional to the magnetic flux, and 

3. the underlying manifold defined by the solution does not include curvature singularity. 
The stationary metric ( |2.2| ) compatible with rotationally symmetric solutions becomes 

2 



ds 2 



dt + T-^rd9 



e 2 v 2 



dr 



K 



e 2 v 2 



b(r)(dr 2 + r 2 d6 2 ), 



(4.1) 



where r = \/x i x l . With the help of a gauge transformation, ansatz for the scalar field is 

4> = vf(r)e~ ine , (4.2) 
and thereby a component of the metric b ( p. 21 ) becomes 

-f2 e -(f 



6(r) 



„2n 



G 



(4.3) 



where F is an undetermined constant due to the harmonic function h(z). Under this (r,6)- 
coordinate, the equations of motion ( |2.23|) - (|2.24j) become 



Id d , f 2 



FG 



r 2n 



G 



(/ 2 -l)[/ 4 -2(l + ^)/ 2 + 0y, (4-4) 
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Id d 

— — r — mip 

r dr dr 



FG 



r 2n 



G 



(/ 4 - 2/ 2 + 



(4.5) 



To obtain nonsingular solution at the origin, Eq. ( |4.2| ) fixes the boundary condition of the 
scalar field at r = 

nf(r = 0) = (4.6) 
which implies the leading behavior of scalar field for small r 



/ ~ /or", 

and Eq. ( |4.5|) gives the boundary condition of ip at r = 

ij)(r = 0) = ip , 

where ipo is an arbitrary positive constant. 

Rewriting the magnetic flux ( |2.11| ) in terms of / and ip 



(4.7) 



(4.8) 



$ = ±- 



2tt 



d\nf 1 

-1 



a — i' — ; — — -(1 — f 2 )r 



dr 



d ha. if) 
dr 



(4.9) 



and the angular momentum (|2.12|) as 
J = 



a 



TTK, 

2^ Jo 



r 



2/i 



(r - 2/ 2 + 



7TK d 

— -r — In w 

e 2 G dr 



(4.10) 
(4.11) 



we notice that, from Eq. ( 4.11| ), spinning objects have long range term of ip for large r, 
specifically ip ~ r~ 2a , where a is an arbitrary number and it also contributes to the magnetic 
flux (|4.9| ) if (poo 7^ 1. The energy (|2.9| ) can also be expressed in terms of / and ip such as 

- roc ( /df\* F /f2 fi -(/a-l) 7/; ^-l^ G 

E 



Kr ~lo drr \Kdr) + 4 



~2n 



G 



/ 2 (/ 2 -l) 2 -^(/ 4 -2/ 2 + 



t) 2 



|e$|. 



(4.12) 



Eqs. ( |4. 12| ) - ( |4.10| ) give possible boundary conditions of the scalar amplitude at r = oo for 
a given a 

|0| _ f or 1, when G[n + a(l - 0^)] < 1 
v arbitrary, when G[n + a(l — 0^)] > 1. 



(4.13) 
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About the geometry of two dimensional space £ described by the 7^ part of Eq. (|2.2|) 
(or by the b(r) of Eq. ( [4.1|) ), various manifolds are characterized by the area of space defined 
by A = 2n / drrb(r). Suppose that the asymptotic behavior of the scalar amplitude as 
\4>\/v ~ r~ £ where e is related to a, the above categories are distinguished by the value of 
G[n + e + (1 — (jx^ja], i.e., whether it is less than or equal or greater than 1. Additional 
quantities characterizing the property of £ are the radial distance from the origin p(r) = 
(e 2 f 2 /|/t|) J*q dr ' Jb(r') and the circumference Z(r) = (2ne 2 v 2 /|/t|)ry / 6(r). When e = 0, the 
above expression of the scalar amplitude \(f>\/v ~ r~ £ corresponds to a nonzero finite (f)^. 
Rapid decay of the scalar amplitude is depicted as e — > 00. So this assumption may not lose 
any generality. In this section we will obtain the rotationally symmetric solutions carrying 
magnetic flux (or charge) and spin. Underlying spatial manifold £ will comprise three types; 
cone, cylinder, and two sphere. We address this question separately for these three cases. 

From now on, let us examine precisely all possible rotationally symmetric solutions of 
the Bogomolnyi equations ( |4.4|) -( fO)|) . When = 1, Eq. (|4.4|) for / is decoupled from the 
equation for ip ( f4.5| ). Thus it may be convenient to investigate the soliton solutions for the 
cases separately by whether (p^ is equal to one or not. 



(a) 



We consider the case 



1 by dividing it into two categories: Gn 7^ 1 and Gn = 1. 



When Gn is smaller than one, we introduce variables u and R such that 



u = \nf 



(4.14) 



and 



R 



,1-Gn 



1 - Gn 



Then the Bogomolnyi equation ([4.4j ) is rewritten as 



(4.15) 



d 2 u 
dR 2 



dU, 



off 



1 du 



du R dR ' 



(4.16) 



where f/ e ff is given by 



U eS (u; (poo = 1) 



F 

— exp 
G 



-G(e u -u-l) (e n - 1 



(4.17) 
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As shown in Fig. ||-(a), U e s\ 



1) has two local minima at u = ln[(G + 1 ± y 2G + 1)/G] 



and three maxima at u — 0, ±oo. If we interpret u as the position of a hypothetical particle 



with unit mass and R as time, Eq. ( }4.16| ) is Newton's equation for 1-dimensional motion 
of the hypothetical particle moving in the potential C/ e g and subject to a time-dependent 
friction, —(du/dR)/R. For n/0 cases, the particle also receives an impact at R = from 
the delta function term, i.e., [2n/(l - Gnf/^-^^R 1 ^ 1 -^) from Eq. ( ggg ). The energy 
of the hypothetical particle £(-R) is defined by 

2 



£(J2) 



1 

2 



(4.18) 



and it is a monotonically decreasing function of i? because of the friction. When n — 0, 
trivial solution k = is the unique solution which describes two dimensional flat space. We 
now show that, only when < Gn < 1, there always exists a finite energy solution whose 
base manifold £ is an asymptotic cone. 

For this we prove that, for a suitably-chosen initial parameter /o, we can obtain a motion 
of the hypothetical particle such that it starts at negative infinity with the initial velocity 
given by Eq. ( |4.7|) and stops at u = at R = oo. 
For r — > 0, the behavior of the scalar field turns out to be 



/ 



/or" 
/or" 



1 + f 6 GeCffr* - ^(1 + G)e*3(l6(2 + G> 



_G ,2(1-G) 

Jo 



16 



Ge^/o 2 ' 



512 



1 + G)e 2G / 4G r 4 + ■ 



f \/,; -G 2 F)/ 4G r 4 + -- 

when n = 1, 
when n > 2. 

(4.19) 

Let us fix an arbitrarily-large number R Q for R in Eq. ( f4.15|) . If we choose /o sufficiently 
small so that / G r 2 << 1 then the terms higher than second-order in Eq. ( |4.19| ) can be 
neglected for R < R . If we substitute this approximated solution ( |4.19| ) into Eq. ( |4.18| ) 
then the energy of the particle £(Rq) in Eq. ( |4.18j ) is estimated by 



^(^o 



n 



2(1 - Gn) 2 Rl 



G f 2G 




1 

Rq 



G 2 n 



R 



2Gn 
1-Gn 



(4.20) 



Let us choose Rq and /o such that the square bracket in the right-hand side of Eq. ( |4.20|) 
negative and the absolute value of this second term is larger than the first term. Furthermore, 
by taking / small enough, additional positive contribution from the negative higher order 
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u = ln[(G + 1 - V2G + 1)/GK 




0oo = 1 
0oo = 

= 71 + 1/2^ 



u 



Z o 

(b) 



Figure 2: (a) Typical shapes of U cS : The solid, dashed, and dotted lines correspond to 
0oo = 1,0, and \Jl + l/lnGv 2 . (b) A typical shape of W e g. 

18 



terms does not change the sign of £(Ro). Then S(Rq) is negative, which means £{Rq) < 
U e s(0). Since the particle energy £ decreases further as R > Ro becomes large due to friction, 
the hypothetical particles has a turning point for sufficiently small /o and stops finally at a 
local minimum of the effective potential u(R = oo) = ln[(G + 1 — \J2G + 1)/G] after some 
oscillation around it. For the large fo and the sufficiently small Rq which makes Taylor 
expansion in Eq. ( [4 . 1 9[ ) valid, it is easy to notice that, for R > Rq, 



u{R) > 2 In / r™ + ^-(R 2 - R 2 ) + R 2 In R /R., 



(4.21) 



where U' h 
value at R 



M — max_ oo<u < (dU e s / du) > 0. The right-hand side of Eq. (|4.21|) has maximum 

1/2 



Ri 



R% + 
u(Rx) > 



An 



U'[l-G(n+M)] 



n 



and then u(Ri) satisfies, for a small R , 



1 - Gn 



\ An(l-Gn) 
In — 1 



A I 



+ 21n/ . 



(4.22) 



Therefore if we choose f sufficiently large as assumed previously, then u(Ri) > 0, i.e., the 
particle went over the hilltop of the potential at u = 0. From these results, continuity now 
guarantees existence of a vortex solution connecting the boundary values, u(0) = — oo and 
u(oo) = 0, for an appropriate f . This completes the proof and a specific example /(= \<fi\/v ) 
is shown in Fig. [| 

Next we discuss that this solution constitutes an asymptotic cone for Gn < 1 case, 
however there exists no solution for Gn > 1 case. When Gn < 1, the radial distance p and 
the circumference I diverge as r goes to infinity. The solutions of Eq. ( |4.4| ) approach their 
boundary values exponentially 



fnl-UKoiil-G^R), 



(4.23) 



where is a constant determined by the proper behavior of the fields near the origin. Thus 
the asymptotic structure of S is a cone with deficit angle 5 = 1i\Gn which is flat. With the 
help of Eqs. fl4.14|) -( |4.15|) , introduction of a new variable such as 



(4.24) 



rewrites the equation ( f4.5| ) as another Newton's equation for a particle whose position is x 



dR? 



2 v 



1)2 _ I^X 

' RdR 



(4.25) 



19 




r 

Figure 3: Plot of rotationally symmetric solution for 0^ = I. Parameters chosen in the 
figures are: n — 1, F — 1, and G — 1/2 for curved spacetime. The solid and dashed lines 
correspond to \4>\/v and J. 
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The external force in the right-hand side of Eq. (|4.25| ) is always negative so that a particle 
starts at \ = In ipo with zero initial velocity and goes to negative infinity as it decelerates to 
zero for large R. Since the external force term in the right-hand side of Eq. ( |4.25| ) decays 
exponentially for large R, the leading term of asymptotic solution for \ is obtained by solving 



the linearized equation ( |4.25|) . Suppose that ip ~ r (a > 0) for large r, these topological 



vortices carry quantized magnetic flux ( f4.9| ) 



$ = ±^ (4.26) 



and angular momentum (|4.11 



2txk 

J = — (4.27) 
e z G 

Note that the spin is arbitrary in curved space which is different from the result of flat 
spacetime, i.e., J oc n 2 in flat spacetime. 

Let us consider the case that Gn is larger than 1. Suppose that there exists a solution 
when Gn > 1. Then the radial distance p(r = oo) is finite and the circumference l(r = oo) 
vanishes since the metric b(r) decreases faster than 1/r 2 for large r, which means that 
the spatial manifold £ is bounded. Since the Euler number given in Eq. ( |2.26| ) must be 
nonnegative due to the fast decaying property of the solution, two dimensional sphere S 2 is 
unique candidate. From now on we call the point which corresponds to r = "the south 
pole" on S 2 and that which corresponds to r = oo "the north pole" on S 2 . The configuration 
of the scalar field <ft & t the north pole does not vanish as <f)(r = oo) = ve me . Therefore the 
scalar field </> is not well-defined at the north pole and there is no regular solution in this 
case. 

Let us study the case when 0^ = 1 and Gn = 1. We introduce R as 

7? = lnr, (-cx)< J R<oo) (4.28) 

which reflects the scale symmetry (r — > Ar) of the Bogomolnyi equation ( f4.4|) at this critical 
value. Then Eq. fl4.4j) takes the same form as Eq. ( [4.1 6|) without friction term. Similar to 
the previous cases, a hypothetical particle moves under a conservative force only and hence 
the Bogomolnyi equation can be integrated to a first-order equation ( |4.18| ) whose energy 8 
of the hypothetical particle is conserved in this case. The initial particle energy S(R = — oo) 
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is determined by the initial behavior of / 

£ 



1 
2 



du 

Jr 



In' 



(4.29) 



H=-oo 

Here if we recall that three maxima at u — 0, ±oo are degenerate, we can easily notice that 
there is no such particle motion which satisfies u = as R — > oo. Therefore it completes the 
nonexistence of soliton solution when ^ = 1 and Gn = 1. 



When 7^ 1, / and ^ are described by the coupled equations ([14] ) and (|4.5| ). As has 
done in the case of ^ = 1, we introduce a new variable 



l-#|JlnV> 



with m and in Eqs. (|4.14j ), ( (4.15 ) and ( [4.28 ). Then the Bogomolnyi equations 
become 



d 2 u 

d 2 tj 
dR 2 



dV e f{ 1 du 
du R dR ' 
d(TVcs) 1 



where =F in Eq. (|4.32| ) denotes — for (p^ < 1 cases and + for (p^ > 1 cases. V e s is 



.£! e -G*(e"-M-l) ( ' e 2n 

G 



2e u + 



(4.30) 



(4.31) 
(4.32) 



(4.33) 



Note that Eq. fl4.31|) contains an impact term at a boundary value of R which corresponds 
to r = if n ^ 0, and the friction term in both equations (|4.31 )- (|4.32 ) vanishes if Gn = 1, 
which reflects the scale symmetry (r — > Ar) of the Bogomolnyi equations (Q)- (fO| ). 

We continue to use the terminology of Newton's equation for the two-dimensional motion 
(u, £) of a particle when (p^ > 1 and for one-dimensional motion of two interacting particles 
u and £, when (p^ < 1. The boundary values of u at R = oo can be read from shapes of the 
effective potential C/ e g. In order to constitute a soliton configuration, minimum requirement 
to the boundary value of the scalar field, v(poc, at spatial infinity is to be an extremum 
of the effective potential U e s. Specifically, from dU e s/du = and e u ( R= °°^ = (p 2 ^, we can 
easily read that possible values of (p^ are 0, 1, \Jl + 2/G and oo. Since (p^ = 1 case was 
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already analyzed, let us discuss the other cases. One point on the cosmological constant 
should be noted: We recall that 0oo = 0, 1, 1/VtcGv 2 make the potential vanish, that is, 
V(\(p\ = v^oo) = 0. If they exist, regular solutions satisfying the other boundary conditions, 



1 + 2/G or 0oo = oo, imply self dual solitons in curved spacetime with a nonvanishing 
cosmological constant. 

(b) 0oo=O 

For this case the scalar field approaches the symmetric local minimum for large r as 
shown in Fig. [TJ (dashed line). Then the Chern-Simons gauge field is topological field, 
i.e., it is no propagating degree, and charged meson is massive, m meS on = (e 2 v 2 /2\n\). All 
the solitonic excitations are the nontopological solitons which are marginally stable since 
E = m meson Q = (e 2 v 2 /2\k\)Q. We consider those solutions with nonzero vorticity {n ^ 0) 
separately from those solutions with no vorticity {n = 0). 

(b-i) n = 

As mentioned previously, we regard Eqs. ( ]4.31| ) - (|4.32|) as Newton's equations for two 
particles of which the positions are u(R) and £,{R), respectively. Here nG = and then R is 
in fact equal to r as shown in Eq. ( |4 . 1 5| ) . Since the force, —dV e s/du, applied to one particle, 
of which the position is u, is negative definite for u in the region — oo < u < and for all £, 
this particle starts out at time zero, R = 0, at a point uq (— oo < u$ < 0) and approaches to 
negative infinity u = — oo, at infinite time R = oo. Similarly, since the force — d(— V e g )/<9£ 
applied to the other particle, of which the position is, £ is positive definite for any £ and u 
in the region — oo < u < 0, this particle initially starts out at a point £o and finally goes to 
positive infinity, £ = oo. 

Near the origin power series solutions give 

fir) « 2o(l-^o (1+6) (l-#o 2 )[2(l + ^^ (4-34) 
tf(r) w ao\l + ^Gg 2 (1+d \2-g 2 )e^ 1 -^a^r 2 + •••!. (4.35) 



Eqs. ( |4.4p -(fP)D do not constrain both g (0 < go < 1) and ao, however if g > 1, both f(r) 
and ip(r) will be monotonically-increasing functions of r and hence the boundary condition 
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at r = oo cannot be met. It is consistent with the argument given in terms of Newton's 
equations for the region of uq (—00 < u < 0). Eqs. ( |4.4j )- (|4~5| ) imply that long distance 



behaviors of the scalar field / and the off-diagonal component of metric if) should be f(r) 



and if>(r) 



„2» 



so that Eq. (|4.3| ) gives b(r) ~ r 2G<y£ a \ Positivity of the energy (|2.10|) 



restricts the values of a and e to satisfy e — a > 0. The precise forms of power series solutions 
for large r are 



f(r) 
if>(r) 



9ooT 



,2a 



1 



1 - 



F(l + G)e 



G 



8[G(e - a) + e - I} 2 00 

FGe G 
[G{e - a) + e - !]'■ 



9 l l+C!) air-^e-a )+ e^ + 



2{1+G) a G r -2{G(e-a)+e-l) , . 

lJ OO OO 



(4.36) 
(4.37) 



where and are the constants determined by the behavior of the fields near the origin. 
It is of our interest to determine the lower bound of G(e — a) + e — 1 by examining the 
correspondence between the short- distance and long-distance behaviors of the solutions, i.e., 
(^O)tto) and (e, a). Since / becomes small near the origin for a sufficiently-small go and an 
arbitrary ao, we can expand the right-hand side of Eqs. ( |4.4j )-( [4T5| ) for the leading order of 
/. If we define g as g(r) = Fe G (l + 2G)f 2( - 1+GS) ip G , then the approximate equation for g is 
the Liouville equation. An exact solution satisfying the boundary condition at the origin is 



g{r) 



go 



2(1+G) G 



a: 







(i + \gT +G) *$r 2 ) 



2\2 



(4.38) 



which is small enough for all r and approaches to the exact solution of Eqs. ( |4.4p - (|4T5|) as 
g — > 0. Comparing the solution ( |4.38| ) with those in Eqs. ( |4.36| )- (|4~3~7l) , we obtain a lower 
bound for e and a, i.e., G(e — a) + e — 1 > 1, and it is consistent with regularity of the power 
series solutions fl4.36j)-([Q7j). As go — > 0, the magnetic flux becomes $ = ±(2 — e)/2eGv 2 



and it implies a constraint for e, i.e., e < 2. 

Furthermore, another category in examining the solutions is to know that the area A 
of spatial manifold £ is finite or infinite. When G(e — a) < 1, the area A, the radial 
distance at infinity p(r = 00), and the circumference at infinity l(r = 00) are infinite. 
Then the asymptotic structure of S is also a cone with deficit angle S = 8tt 2 Gv 2 (e — a). 
Though the global structure of the cone is the same as that of vortices of (p^ = 1 case, these 
manifolds have long-range tails for large r as shown in Fig. |§-(a) and the shape of the angular 
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momentum density (see Fig. §-(b)) shows that it is accumulated near the origin. Since the 
circumference at infinity l(r = oo) is finite when G(e — a) = 1, the asymptotic structure 
of £ is a cylinder. When G(e — a) > 1, the space S is bounded manifold and then a two 
sphere is unique candidate. The Euler number fixes G(s — a) = 2 that is b(r) ~ r -4 . Then 
the regularity at the north pole (r = oo) requires that the scalar field behaves as / ~ s n for 
s ~ 0, where s is radial coordinate in the new coordinate system whose origin is at the north 
pole. However, since the scalar field <fr does not vanish at the south pole which corresponds to 
the point r = or equivalently s = oo, it is not well-defined. Specifically, <ft — \4>\e me = g§e 
at the north pole. Therefore there is no regular solution when G(e — a) > 1. 

Since the solutions with |0|(oo) = are nontopological solitons, they are characterized by 
the U(l) charge (or equivalently the magnetic flux) Q{= — k$) = =(=(27r«;/e) e+ (e 2 G/2iiK) J 
and the angular momentum J = —(2irK/e 2 G)a, which need not be quantized. Note that the 
sign of angular momentum for these nontopological solitons is opposite to that of vortices 
as given in Eq. (|4.27|) as in the case of flat spacetime. In case of the cylinder the magnetic 
flux has a fixed value $ = 27rn/eG for a given set of parameters. 

(b-ii) n ^ 

For n / solutions, the motion of a hypothetical particle whose position is depicted 
by £ coordinate resembles that of the n = solution since it follows the same equation of 
motion ( |4.32| ) and the starting point is also an arbitrary constant £ - O n the other hand, 



the motion of the other hypothetical particle whose position is u coordinate is different from 
n = solution due to the vorticity. Specifically, the particle should start out at negative 
infinity, u = — oo, with an initial velocity given by Eq. ( |4.7| ) which is induced by receiving an 
initial impact at time corresponding to r = 0, and then it turns at an appropriate position 
"Umax at a certain time, and finally returns to u — —oo at time R(r)\ r= o. We now argue that 
there exist such nontopological vortex solutions by showing the existence of turning point 
w max between — oo and 0. 

For small r, power series solutions are 

f(r) « V" [l- T/ ^ = V2 (1 + G)e^f +5) /? V(" +1 > + ■■■], (4.39) 
L 8(n + iy J 

V>(r) w foil + FGe d h 2 {1+d) $ r 2 ( n+1) + ■■■}. (4.40) 
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Figure 4: Plot of rotationally symmetric nontopological solitons for 0^ = and n — 0. 
Parameters chosen in the figures are e — a — 1, F = 1, G = 1/2 for curved spacetime: 
(a) \<f>\/v vs f (b) J vs f . The solid and d2§hed lines correspond to /(0) = 0.8 and 0.3, 
respectively. 



Note that the sign of second-order terms is negative for / and positive for ip, which looks 
consistent with the above description in terms of Newtonian mechanics. Let us assume 
an arbitrarily-large number R for R both in Eq. Q4.39D and Eq. ([4.40 ). If we choose ho 
sufficiently small for a given (3q, then the terms higher than second-order both in Eq. (|4.39|) 
and Eq. Q4.40Q can be neglected for R < R Q . 

The leading terms of / and ip f° r large r are assumed to be f(r) ~ r~ e (e > 0) and 



2a 



a > 0) so that we obtain a condition n + e — a > from positivity of the energy. 



Long distance behavior of / and ip is 

F(l + G)e d 



f(r) 



hr^r' 



1 



2a 



1 - 



8[G(n + e - a) + e 
FGe d 



1] 



_ h 2(l+G) iG -2(G( 



n+e— a)+e— 1) 



+ ■ ■ 



_j 1 2(l+G) pG r -2(G(n+e-a)+e-l) + 



4.41) 

, (4-42) 



[G(n + e-a) + e- l] 2 x 

where h^ and (3^ are the constants fixed by requiring the proper behavior of the fields near 
the origin. Through an approximation to the Liouville equation as was done in n = case, 
we obtain a condition G(n + e — ct) + e — 1 > n + 1 for each n and the equality holds when 
h$ — ► 0. In this limit the magnetic flux become $ = ±(n — e + l)/2eGv 2 which implies the 
upper bound of e, i.e., e < n + 2. 

Analysis of the asymptotic structure of spatial manifold E is almost the same as n = 
case, if we replace e— a to n+e— a. The n^O solitonic configurations support asymptotically 
a cone with deficit angle 5 = 87r 2 Gv 2 (n + e — a) when G(n + e — a) < 1 and a cylinder 
when G(n + e — a) = 1. On these open spaces, the nontopological vortices are hybrids of 
the vortices at short- distance region and the nontopological solitons at long-distance region 
as shown in Fig. 5-(a)(b). 

When G(n+e — a) > 1, the space £ is bounded and then it should form a two dimensional 



sphere. Then the Euler invariant (|2.26| ) should also be equal to that of a smooth S 2 such as 

G(n + e-a) = 2. (4.43) 

Now we show that both e and a can not be arbitrary by imposing the regularity condition. 
From the behavior of the metric b(r) at large r, the radial distance p(r) behaves near the 
north pole as p — p(oo) ~ r -G{n+e-a)+i _ Q n other hand, if we choose a new radial 
coordinate t around the north pole as p — p(oo) ~ t, then the regularity demands that the 



27 




200 400 600 800 

f 

(a) 




100 200 300 400 500 600 

r 



(b) 

Figure 5: Plot of rotationally symmetric solutions for 0^ = and n — 1. Parameters chosen 
in the figures are e — a = 1, F — 1, G — 1/2 for curved spacetime. (a) \4>\/v vs f (b) 
J vs f. The solid and dashed lines correspond to (/loc^oo) = (0.02,0.2) and (0.0017,0.2), 
respectively. 



scalar field behaves as / ~ t n for small t. Comparing these with Eq. ( |4.7|) , we get e = n 
and a = 2(n — Moreover, since ip expressed in t coordinate should also be a positive 
constant at the north pole, a (equivalently the angular momentum J) should vanish and 
thereby the solution can exist only when n = 1/G. For this sphere case with Gn = 1 the 
Bogomolnyi equations in Eq. ( |4.31 )-( 4.32 ) have no friction term. Therefore, one of those can 



be integrated to a first-order equation 



l(^Y-im\ v ^ M = 2A (4.44) 



2\dRJ 2\dRJ 

where the integration constant in the right-hand side is determined by considering the be- 
haviors of u and £ near r = in Eqs. (|4.39| )- (|4.4CI| ). This sphere is symmetric under the 



inversion between the south pole and the north pole so that it can also be interpreted as a 

configuration that two vortices with vorticity n lie both at the south pole and at the north 

4vr 1 

pole. The total magnetic flux is twice of the flux of a hemisphere, $ = ± — n = ±— — -. 

e Gev 2 



( c ) <Poo = yi + i 

Suppose that for each n there exists a solution to satisfy this boundary condition. Since 
the shape of effective potential U e s given by the dotted line in Fig. 0-(a) may guarantee the 
existence of a solution, the finite energy condition in Eq. Q4.12 ) restricts G(n+a(l— 0^)) > 1 



and it implies the volume of two dimensional space £ should be finite. Two sphere is the 
unique candidate as explained in 0^ = case. Since 0oo ^ contradicts to the condition 
for the regularity on the north pole of S 2 , there does not exist such regular solution. 



(d) 000 = 00 

If we look at the scalar potential unbounded below as given in Fig. [l] and an expression 
of the energy in Eq. ( f4.12| ), we notice that there is no reason to exclude (j)^ = oo solution as 
a candidate of a finite-energy solution. When n = 0, u = oo is the solution which describes 
the potential minimum at negative infinity. However the energy of it is of course negative 
infinite so that it is the solution out of our scope. When n^O, the conservative force by the 
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effective potential, dU e g / du, is singular at u = and thereby there is no regular solution to 
connect both boundaries u = — oo and u = oo. 

5 Behaviors of Multi-soliton Solutions 

In the previous section we restricted our interest to self-dual solitons with rotational sym- 
metry. We study the behaviors of arbitrary multi-soliton configurations without rotational 
symmetry in what follows. 

First, we show the fast decaying property of topological vortices on asymptotic region 
irrespective of both parameters of the theory and positions of vortex centers. Then nonex- 
istence of the topological multi- vortices are proved when Gn > 1. 

We prove the fast decaying property of solutions of Eq. ( |2.23| ), which is not necessarily 
radial. 

Theorem 1: Assume that a solution f 2 for Eq. ( |2.23| ) satisfies the following conditions: 
< f 2 < 1, lim f 2 = = 1, < Gn < 1, and e h+h = cq is a positive constant. Then 

2— >00 

f 2 converges to 1 faster than any polynomial rate as \z\ increase in the following sense: f 2 
satisfies \f 2 (z) — 1| < ■At as \z\ — > oo for any positive k. 

Next we show the following nonexistence of multi-vortices. 

Theorem 2: Assume that < f 2 < 1, Gn > 1, lini^oc f 2 = 0^ = 1, and e h+h = 
Co is a positive constant in Eq. ( |2.23| ). Then there is no solution f 2 of Eq. ( |2.23| ) with 
/dz|/ 2 -l| 2 <oo. 

Detailed proofs of Theorem 1 and 2 are presented in Appendix. 

6 Conclusion 

In this paper we have precisely described the self-dual Chern-Simons solitons in curved 
spacetime. Bogomolnyi type bound for the original and the dual-transformed theory has been 
derived under a one-parameter family of 8 scalar potential. In the context of the duality 
transformation in continuum, the role of Higgs and the topological sector of scalar phase have 
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been clearly demonstrated by showing how they support the nonperturbative excitations, 
specifically the soliton spectra, and affect to their mutual interactions through introduction 
of a Jacobian in path integral measures. Then we analyzed the Bogomolnyi equations under 
a rotationally symmetric ansatz and obtained all possible soliton solutions. Although the 
scalar potential has an additional </> 8 term with negative coefficient in comparison with that 
in flat spacetime and changes its form along the boundary value of the scalar field, this 
system supports regular soliton configurations with finite, positive energy whose boundary 
values are the same as those in flat spacetime, i.e., |0| (oo) = v or |0|(oo) = 0. In the 
former case there are topological vortices which carry a quantized magnetic flux (or a U(l) 
charge) and constitute an asymptotic cone and cylinder as the underlying spatial manifold. 
The latter case contains the nontopological solitons and vortices which have a continuous 
variable charge and whose underlying manifolds comprise an asymptotic cone, an asymptotic 
cylinder, and a two sphere as the value of total magnetic flux is less than or equal to a critical 
value. All of them are spinning objects whose values can be arbitrary in contrary to the 
vortices in flat spacetime. 

While the existence and the properties of rotationally-symmetric self-dual solitons have 
been rigorously examined, the stability of those solutions is an open question. Perhaps the 
unboundedness of scalar potential can be understood by a supergravity version of the model 
as the supersymmetry did in flat spacetime pi] , but the stability of nontopological objects 
outside the Bogomolnyi limit both in flat and curved spacetime may await investigation 
using numerical computation. In relation with the cosmological constant, asymptotic region 
of every open space formed by the self-dual solitons has zero cosmological constant though 
there exist several local vacua with nonvanishing cosmological constant. It may imply a 
connection between the self-dual Chern-Simons solitons and supergravity theory 123]. This 
seems also consistent with unattainability of a Bogomolnyi type bound in a gauged system 
with nonvanishing cosmological constant [p3fl . For multi-soliton solutions, which are not 
necessary rotationally symmetric, fast decaying property is proved and nonexistence of finite 
energy solutions are presented when Gn > 1. Though some aspects of solitons in dual- 
transformed theory were briefly commented, the quantum field theoretic issues such as the 
phase transition induced by Chern-Simons solitons need further study. 
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7 Appendix 



In Appendix, detailed mathematical proofs of Theorem 1 and 2 are provided. We denote 
dz = d Xl dx 2 , A = d 2 , a = i+g-V^ , d = 1±g±^5 and 



where 



and 



H(z,T) = H 1 (z,T)+H 2 {z,T), 



Hi(z,T) = A In T 



Gen T G e G ( 1 " T ) 
F 2 (* T) = — — -(T — 1)(T — a)(T — d). 

When e* 4 * = c and 0^ = 1, Eq. (ggg ) is equivalent to F(z,T) = 0. Let H 2 (R 2 ) be 
the Sobolev space, which is the completion of C^°(R 2 ) with respect to the norm \\w\\ = 
{J R 2 \Vw\ 2 + w 2 dz)*. 

Proof of Theorem 1. 

Take u(z) = 1 — c/\z\ k . We will show that u < f 2 on the outside of a large ball. 
Step 1. H(z, u) > on the outside of a large ball. 
Take a fixed constant b so that a < b < 1 and 

^H 2 (z,T)<0, (A.l) 

when 6 < T < 1 and |z| is large. Since f 2 goes to 1 at infinity, we can take large ro so that 
f 2 (z) > b on \z\ > r and we take c so that w(r ) = 1 — c/tq = 6. Note that 

ck 2 \z\ k ~ 2 



A In w 



(\z\ k -c) 2 

ck 2 
\z\ k+2 u 2 



We calculate H(z, u) on \z\ > ro- Since \u — a\ > \b — a\ and \u — d\ > \d — 1|, 

C (5 M G e (l-«)G 

#2(2, m) = - /rrP=n . — T ,(u- l)(u-a)(u-d) 



2 (n^ik-^i 2 ) G ' 
- #MS^W (6 " a,(<i " 1) - 
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Let us denote a = G{b — a)(d— 1)cq/2. Since (Ilp=i \ z 
d and \z\ > ro, we estimates H(z,u) 



\ 2 ) G < c'\z\ 2Gn 



for some constant 



H(z,u) > 



> 



c 



\k+2 
C 



au 



G 



> 



I fc+2 
C 



a 6' 







2 






2 


2\G 


(n£=? I 






2-jG 



- fc 2 /6 2> ) 



I fc+2 



(a^lz^ 1 "^ - k 2 /b 2 ) . 



(A.2) 



Eq. ( |A.2|) is nonnegative when \z\ > r$ by taking b which is sufficiently to 1. 

Step 2. f 2 (z) > u(z) on D = R 2 — B(r ), where r is defined in step 1. 

Assume that f 2 = e w and u = e m on D, where w and m are smooth functions on R 2 . 
Note that H(z, f 2 ) = H{z, e w ) = 0, H(z, u) = H(z, e m ) > on D, and f 2 (z) > u(z) = b on 
dD. From Eq. (|A.2|) , we have 



> H(zJ 2 )-H(z,u) 
= H(z,e w ) - H(z,e m ) 
= A(w-m) + H 2 (zJ 2 )-H 2 (z,u). 



(A.3) 



From Eq. ([ATI 



^(^TO - ff 2 (*,r 2 ) = -A(z)(Ti - T 2 ) 
for some nonnegative function X(z) for 6 < T 1 ,T 2 < 1. Therefore, 

A( w - m ) < -# 2 (z, / 2 ) + H 2 (z, u) = X(z)(f 2 - u). 



(A.4) 



Assume that f 2 — u < on a domain £)' C -D, then A(u> — m) < on D'. By the Maximum 
Principle w — m cannot have minimum inside of D' . Since w = m on dD' and w < m on 
D £>' must be the empty set. We conclude that f 2 > u on D. We proved Theorem 1. 

Assume that there exists w G Hf(R 2 ) so that the solution f 2 = e w < 1 for Eq. (2.22) on 
D. We claim that f 2 converges to 1 as \z\ — ► oo. 
From Eq. (|2~23| ), 



Aw = -H 2 (z,e w ). 
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For a point q G D, \w(q)\ < C (/b((j,i) \ w \ 2 ^{z, e w ) 2 dz*j by the standard elliptic estimates 
for w on R 2 (see [[24]]). Since < T = f 2 = e w < 1, J B{q l) H 2 {z, e w ) 2 dz goes to as \q\ -> cx). 
Since w G Hf(R 2 ), \w(q)\ — > as |g| — > 00. Therefore, / 2 (-z) ^ 1 as \z\ — ► 00. From the 
above we proved the following Theorem. 



Theorem A: Assume that Gn < 1, 0oo = 1, e h+/l = c is a positive constant, and 



there exists w G H 2 (R 2 ) so that a solution f 2 = e w < 1 for Eq. (ggg) on D = R 2 - B^), 
where D is the outside of some large ball of radius R\. Then f 2 converges to 1 faster than 
any polynomial rate as \z\ increases in the following sense: f 2 satisfies \f 2 {z) — 1| < as 
\z\ — >• 00 for any positive k. 

Proof of Theorem 2. 

We use the same notations as in the Proof of Theorem 1. Take u(z), b, and ro as in the 
Theorem 1. Note that f 2 {z) > b when \z\ > r . Choose b' so that b < b' < 1. Take r\ > r 
so that f 2 (z) < b' on \z\ = ri. Take c so that u(r\) — 1 — cjr\ — V . Since \u — a\ < |1 — a\ 
and \u — d\ < \d — b\, 



H 2 (z,u) 



c G ifie^-^ 6 
~ (Tl^i \z-z p \ 2 f 



(u — l)(u — a){u — d) 



< 



CqGc 



u 



G 



2N fc (nPk-^| 2 ) G 



-(1 -a)(d-b). 



Let us denote a' = G(l — a)(d — b)c /2. Since (Ilp=i \ z ~~ z p\ 2 ) G > c\z\ 2Gn for some constant 
c and \z\ > r±, we estimate H(z,u), 

I -vl2 



H(z,u) < 



< 



< 



lfc+2 
C 



|fc+2 
C 



a'u G 



a 





2 


2 











5 - k 2 /u 2 

k 2 



\k+2 



(a'c\z\ 2 ^- Gn ^ - k 2 ) 



(A5) 



Since Gn > 1, Eq. (|A.5 ) is nonpositive when \z\ > r\ by taking b' which is sufficiently close 
to 1. 

Next we show that f 2 (z) < u(z) on D = R 2 — B(ri), where r± is defined in above. 
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Assume that f 2 = e w and u = e m on D, where w and m are smooth functions on R 2 . Note 
that H(zJ 2 ) = H(z,e w ) = 0, H(z,u) = H(z,e m ) < on D, and f 2 (z) < u(z) = b' on 3D. 
From above, we have 

< H(z,f 2 )-H(z,u) 

= H(z,e w ) - H(z,e m ) (A.6) 
= A(w-m) + H 2 (z,f 2 )-H 2 (z,u). 

From Eq. (TAD , 

H 2 (z,T 1 ) - H 2 (z,T 2 ) = —X(z)(T 1 - T 2 ) 

for some nonnegative function X(z) for b < Ti,T 2 < 1. Therefore, 

A(w-m) > -H 2 (zJ 2 ) + H 2 (z,u)=\(z)(f 2 -u) (A.7) 

Assume that / 2 -ti>0ona domain D' C D, then A(w — m) > on D' . By the Maximum 
Principle w — m can not have maximum inside of D' . Since w = m on 3D' and w>mon 
D', D' must be the empty set. We conclude that f 2 <u=l — c/\z\ k on D. 

For w G H 2 (R 2 ), J R2 \e w — l\ 2 dz < cexp(ci J R \Vw\ 2 + w 2 dx) for some constant c and c\ 
(see HU). Therefore, \e w — 1| 2 = \ f 2 — 1| > c 2 /\z\ 2k cannot hold for k < 2 as |z| — > oo. Thus 



we proved Theorem 2. 

Remark When (p^ = 1, A In ip < by Eq. (2.23). Since there is no lower bounded 
super harmonic function on R 2 , there is no constant c such that ip > c 2 > and ip has no 
local minimum. 

Next we study the behavior of f 2 when 0^ = and < f 2 < 1. When (p^ = 0, 
A In if) > by Eq. ( |2.24| ). Since there is no upper bounded subharmonic function on R 2 , 



there is no constant c such that ip < c 2 and i/j has no local maximum. From Eq. ( [2.24| ), we 
obtain a decaying condition for f 2 in the following. 



Claim Let f 2 be a solution for Eqs. (|2.23|) -( [2~2"4D when e h+h = c is a positive constant 



2(nG-l) 

and 0oo = 0. There is no positive constant c such that / > c\z\ ^ s for large \z\. 
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Proof of Claim: Take ip G = e h in Eq. Q2.24 ). For large \z\, Eq. (|2.24j) turns into the 
following equation 

n2 f 2G G(l-f) h 

= G c / e ^_ /2(2 _ /2) (A g) 

f2+2G e h 

£ to-^ (A - 9) 

/•2+2G 

where c and c' are some positive constants. Since there is no solution of the equation 
Aw = K(z)e w when K(z) is positive and K(z) > c'\z\~ 2 as \z\ — > oo (see [25[]), we conclude 



2(nG-l) 



that there is no solution for ^ when / 2 > c|z| !+ G for any positive constant c and large 
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